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In this paper, we present an invariant perturbation theory of the adiabatic process based on the concepts 
of t/(l)-invariant adiabatic orbit and (/(l)-invariant adiabatic expansion. As its application, we propose and 
discuss new adiabatic approximation conditions. 

PACS numbers: 03.65.Ca, 03.65.Ta, 03.65.Vf 



Since the establishment of the quantum adiabatic theorem 
[1,2,3,4] in 1923, many fundamental results have been ob- 
tained, such as Landau-Zener transition [5], the Gell-Mann- 
Low theorem [6], Berry phase [7] and holonomy [8]. Also 
the adiabatic processes find their applications in the quantum 
control and quantum computation [9,10,11,12]. Recently the 
common-used quantitative adiabatic condition [15,16,17] has 
been found not able to guarantee the validity of the adiabatic 
approximation [13,14]. Consequently various new conditions 
are conjectured and a series of confusions and debates arise. 
For example, it was argued [18] that the traditional adiabatic 
condition did not have any problem at all and that the in- 
validation of the condition did not mean the invalidation of 
adiabatic theorem [19]. Some new conditions proposed in 
[20,21] but too rigorous to be used conveniently. Although 
[22] also adopted the adiabatic perturbation expansion but did 
not give out proper condition because the basis in [22] can 
not show certain geometric properties in the adiabatic process. 
[23] pointed out the limitation of traditional condition but also 
did not give out a proper condition. To solve the problem 
of insufficiency of traditional adiabatic condition in [13,14] 
and clarify the subsequent confusions, we introduce the con- 
cepts of adiabatic orbit, U{l)-invariant adiabatic orbit and 
U{V)-invariant adiabatic evolution orbit. The meanings of 
adiabatic evolution is reclaimed. And new adiabatic approx- 
imation conditions based on the U{Y)-invariant adiabatic 
expansion with the time-dependent coefficient are proposed 
and illustrated by two examples. 

Let us consider a quantum system governed by a time de- 
pendent Hamiltonian Hit) and the initial state of the sys- 
tem is an eigenstate |m, 0) of H(Q) with eigenvalue E,„{Q), 
where m denotes the initial value of dimensionless quantum 
number set. By introducing a dimensionless time parame- 
ter T - E,n (0) tjti and a dimensionless Hamiltonian hij) - 
H{t)/E„i (0), the time dependent Schrodinger equation reads 



/ ' " " = h(T)\^„lT)), \^„M) = |m, 0) . (1) 
or 

The exact solution |<1>,„ (r)) to Eq.([T]i is referred to as the sys- 
tem's dynamic evolution orbit in the Hilbert space. 

Furthermore, by considering t as a fixed parameter, we can 
always solve the following quasi-stationary equation of the 



Hamiltonian h (j) 

/!(T)|^„(r)> = e„(T)|^„(T)>. (2) 

And the eigenstate |^„(t)) with the corresponding initial state 
|n, 0) is referred to as the adiabatic solution or the adiabatic 
orbit of the system. 

For convenience, we denote y,™, = / (i^„(r) | (pmij)) and the 
dot here and below expresses the derivative with respect to 
time. Apparently, an adiabatic orbit multiplied by an arbitrary 
time-dependent phase factor still describes the same adiabatic 
orbit. It is not difficult to see that the following adiabatic orbit 

lof '-(t)) = exp |-/ j| \e„AX) - j,n,Am^ l'^m(T)) (3) 
is invariant [26] under the following IJ{\) transformation 

\^nAT)) ^ e'A'W |.^„,(t)) a,(0) = 0). (4) 

Here /„(0) = is because of given initial state. We call 
this adiabatic orbit with special choice of the time-dependent 
phase factor as the U(l)-invariant adiabatic orbit. 

It is clear that, although the initial conditions \in, 0) are the 
same, the dynamic evolution orbit |<1),„(t)) do not always coin- 
cide with the adiabatic orbit |^,„(t)), or they are not even close 
to each other. Obviously they coincide if and only if 

r,„„ = (Vn ^£ in). (5) 

In this case, Eq.(l3]i being the solution to both Eq.([T]i and Eq.(|2| 
is referred to as the U{l)-invariant adiabatic evolution orbit, 
describing a strict adiabatic evolution orbit of the system. 

Generally speaking, the dynamic evolution orbit \^,„ (t)) 
starting from the initial state \m, 0) will change among some 
adiabatic orbits which will cause transitions between different 
them. Our task is to find the proper condition under which the 
dynamic orbit is sufficiently close to the adiabatic orbit when 
the Eq.(|5l) is not satisfied. 

Since the Hamiltonian /i(t) is Hermitian, all the t/(l)- 
invariant adiabatic orbits in Eq.© at a given time constitute 
a complete orthonormal basis of the system. In this basis, the 
dynamic evolution orbit of system reads 
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The expansion in Eq.® is referred to as the U(l)-invariant 
adiabatic expansion with the time-dependent coefficients. 
Therefore, the set of coefficients equations reads 



(7) 



where the diagonal elements of the marix M(t) are zero and 
the non-diagonal elements of M(r) read 

M(t),„„ = / (cD«f' (r) I <hf'(T)) = |y„„(T)| (8) 

where 

Omnir) = I dA ie„(A) - e„{A) + y„„ - rmm)+argy„„(T). (9) 
Jo 

Thus, the probability of staying in adiabatic orbit |'l'°f (t)^ 



Pn, (t) = |c„ (t)P 



where t is time ordered operator And one can obtain further 
detailed analysis on Eq.lfTOli in our another paper [24]. 

Accordingly, the adiabatic approximation of system re- 
quires 
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/ r dAM(A) 
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, (10) 



P„,(t)^1. 



(11) 



It means the transition probability from dynamic evolution 
orbit to other adiabatic orbits (except the adiabatic orbit 
r))) can be neglected. 
According to the perturbation theory for Eq.lfTOli. the first- 
order approximation of Pm(r) is 



Pm(T) = 1 - 



je„M) 



dA 



(12) 



Therefore, the adiabatic approximation requires 



I \ynm(M < 

Jo 



dA 



(V« ^ m). 



(13) 



For general situation \0„,„(t) - 0„,„(O)\ > In and |0,„„(t)| > 1, 
the integral of Eq.(fT3b will be sufficiently small, if the phase of 
the integrated function vibrates fast enough and the amplitude 
of the integrated function is small enough, thus we will have 
following adiabatic condition 



^nm(T) » lr«m(T)l (V« + m) 



(14) 



that is 



k„(T) - e„,(T) -1- A„„(t)| » |r„„,(T)| (yni^m) (15) 



where 

^mn ('^) = ymm (T) ~ 



Inn ( t) + — 



d 

— argy„„ (t) (Vn ^ m). (16) 
dT 



Here A,„„ referred to as quantum geometric potential is a new 
quantity comparing to the traditional adiabatic condition. And 
it should be noticed [26] that the quantum geometric potential 
is also f/(l)-invariant under the transformation Eq.©. It 
should point out that A„,„ appears naturally in our theory. And 
for quantum geometric potential one can obtain further de- 
tailed analysis and application in our another paper [25]. Fur- 
thermore, from Eq.(fT3b. the change of the phase in the in- 
tegrated function should be much larger than the integral of 
amplitude, we will have another condition in integral form 
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dA {e„iA) - eJA) + 



f 
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» I dA |y„,„ (^)l, V n m. 

(17) 

Based on Eq.(|7]) and Eq.®, we also prove [27] a theorem 
related to condition Eq. (fT4l i and Eq.lfTT]). Following analysis 
and our subsequent works [24,25] indicate that Eq.(9-18) are 
of abundant content. 

Next we will give two examples to show the validity of 
Eq.O and Eq.^. 

The first example is to indicate that the problems shown 
in [13,14] do not exist because the relation between system a 
and b stated in Ref.[13,14] does not guarantee our condition 
Eq.lfTSll. The proof is given below. 

[14] showed that for an arbitrary time-dependent sys- 
tem a with Hamiltonian /z°(t) and quasi-stationary equation 
/i''(t) In^ir)) = e"(T) \n"(T)), one can construct time-dependent 
system b with Hamiltonian /z''(t) and quasi-stationary equa- 
tion /!*(t) |«*(t)) = e*(T) |«*(t)) as follows 



h''(T) = iU(T)UHT),h''(T) = itlHrWiT) 
yiT)) = UHt) |n''(T)> , e'jT) = -e';,(T). 

Simple calculation yields 

rin (t) = -e"jT)6„„ + rL (t) • 

Using Eq.lfTSll. we have 



(18) 



(19) 



\rnm(r)\ 



« 1, 'in + m 



\eUT)-e%{T)-^lAr)\ 
for system a, while for system b, we will have 



(20) 



k^(T)-e!;(r)-AL(T) |AL(t)I 



(21) 



Comparing adiabatic conditions between system a Eq.(l20l) 
and system b Eq.(l2TI). the denominator of Eq.(l20l) has an extra 
term e°,(T) -e"(r). Therefore, we can not conclude that system 
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b still satisfies Eq.lfTsTl. This indicates that the problems stated 
in [13,14] will not exist if we adopt the new adiabatic condi- 
tion Eq.(fT5l). And it is worthwhile to point out that we can 
present accurate adiabatic condition for the examples given in 
[14] according to our condition Eq.lfTSll. 

The second example is to consider the well-known model, a 
spin-half particle in a general magnetic field. The Hamiltonian 
of the system is 



h(T) - rjo-. + ^(r) {cTx cos IrfT + cr,, sin 277t) 



(22) 



where rj - fujJo/E±(0), ^(r) = ti(L)(T)/E±{0) and wo is a con- 
stant. Suppose that the initial state of the system is |+, 0) - 
exp(-/cry0(O)/2) l+gj.^ with energy eigenvalues £±(0) at ini- 
tial time respectively, where 0(0) - arctan(w(0)/a;o) and |±ec^ 
are the eigenstates of ctj. The dynamic evolution orbits of the 
system are 



|0±(r)) = e-''^^''^e-''^'£^^^^''^ |+,0) . 



(23) 



Here, the two f/(l)-invariant adiabatic orbits passing through 
the corresponding initial eigenstates |+, 0) are 

jjpacfifl ^^^^ ^ g-inT^:g-iBiT)aJ2^-i<^: £ ^ (24) 

here Q(t) = ^J^iT)^ + rf, cos0(t) = ///^(t) and \(t)f"(0)) = 
l±,0). 

Now, we use the new adiabatic condition Eq.lfTTb to ex- 
amine under what circumstances the evolution of the system 
keeps in the adiabatic orbit, jO"''"' (t)^ First of all, it is easy 
to calculate the probability of finding the dynamic evolution 
orbit of the system in the adiabatic orbit (r)^ 



p„,^\{<:>f''(T)\o,(T))\ =- + - 



2 1 1 ^(0)f (t) -I- 1]^ COS 26 



Q(0)Q(t) 



(25) 



where S=-£ ^(A)dA. 

Suppose that |0-+(t)| > 1 and |6»_+(t) - 6i_+(0)| > In, then 
according to the condition Eq.lfTTll. we will have 



^ Jo " n(A) y 4n(Af-. 



+ 1. 



(26) 



After simple analysis, we will obtain following sufficient con- 
dition 



^(t) » T] 



or 



?7 » ^(r) and I ^(A)dA <sc 1 
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(27) 



(28) 



It is easy to see that when Eq.dZTli or Eq, 
Eq.(l25]l. we have 



is satisfied, from 



P,„ = |((Df«(T)|<D+(r))|'=« 1. 



(29) 



Namely, our new adiabatic condition Eq.lfTTb guarantees the 
evolution of the system is an adiabatic evolution. What is 
more, we can choose proper ^(r) to be periodic, so the condi- 
tion stated in Ref.[20] can not be applied to our example and 
has obvious limitation. In fact, the adiabatic process may be 
a longtime vibration process, so the condition influenced by 
the times of vibration stated in Ref.[20] is too much strict. As 
for the general sufficient condition stated in Ref.[21], it is not 
only too complicated to operate but also too rigorous to apply. 

In conclusion, according to the concepts of adiabatic orbit, 
t/(l)-invariant adiabatic orbit and adiabatic evolution orbit 
stated in our paper, we reclaim the meanings of adiabatic evo- 
lution and present an invariant perturbation theory of adiabatic 
process based on time-dependent L'(l) invariant adiabatic ex- 
pansion. Of course, Eq.lfTSll and Eq.lfTTTi can not be proved 
to be sufficient condition, thus we give out sufficient condi- 
tions in [27] which contains Eq.lfTsl) and Eq.lfTTli. however, 
the second condition in [27] are too strict to exclude many 
interesting physical systems. As far as we know, the condi- 
tions Eq.dTSll and Eq.lfTTli can be used to determine whether 
the evolution of the system is adiabatic or not for all familiar 
examples listed in the existed papers. We also preliminarily 
show the influence of quantum geometric potential in the new 
adiabatic condition. Further detailed discussions on quantum 
geometric potential and the new adiabatic conditions can be 
seen in [24,25]. 
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[26] Here we will prove that adiabatic basis (t)^| and 

quantum geometric potential A,„„ (t) are both invariant under 
arbitrary time-dependent t/(l) transformation. 
Proof: Perform an arbitrary time-dependent (7(1) transforma- 
tion on adiabatic orbit \(p„ (t)) 

Wr)) \v'„(T)) = e'^"*^' \ip„{r)) . (30) 
Thus, we will have 

,-.7;[.»(..)-r;,„„(^)]<'^|^;_(^)) 

= e'/».(0)g-'r[^»,W)-r™W]<'.i|^„_(T-)) ^ ' 

The given phase of initial state |m, 0) requires /j,, (0) = 0. Thus 
we prove adiabatic basis jlo^f'" (t)^| is f/(l)-invariant. And the 
proof for A,„„ (t) is similar. Performing the same f/(l) transfor- 
mation, we will have 

= [ym,n{T) - /„(t)] - [7„„(t) - /„(t)] 

+ f arg [i <v„(T)| K^c,„(t)|>] (32) 

= A„„,(t) (ra ^ n) 



These two proofs show that even if the phase yj-j. (r) is de- 
pendent of time, we still can use them to construct meaning- 
ful quantity which has physical meanings and will not change 
under any time-dependent t/(l) transformation. The quantum 
geometric potential A,„„ (r) and adiabatic basis {lOIlf'" (t)^| are 
just the cases. And all these quantities should exist in all time- 
dependent process. We should note that the traditional Berry 
phase can be included if we perform cyclic integral on the phase 
in{K"'(T))}. 

[27] Theorem For an N -level quantum system and an arbitrary real 
B and a time period T, if the following conditions hold 

max,„„„||fl|<3| (33) 

max„,./;dr |f(^)|<3| (34) 
max„^m_ 

X dT \D,„iT)\ < I (35) 

then the probability of finding dynamical orbit in the adiabatic 
orbit |0°f(T)^ is greater than (1 - ef. 
Proof: From Eq.Q and Eq.l[8]l, we have 



CniT) - 1 = / V r dr \yUr)\e-"''~''^c„{r) 




Then from Eq.(33), Eq.(34) and Eq.(35), we have 

l-|c„,(T)|<|l-c,„(T)|<e (37) 

namely, 

P^ = lc„(T)p > (1 - sf . (38) 
Thus we prove the theorem. 



